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Recent experiments with increasingly larger numbers of qubits have sparked renewed interest in adiabatic 
quantum computation, and in particular quantum annealing. A central question that is repeatedly asked is 
whether quantum features of the evolution can survive over the long time-scales used for quantum annealing 
relative to standard measures of the decoherence time. We reconsider the role of decoherence in adiabatic quan¬ 
tum computation and quantum annealing using the adiabatic quantum master equation formalism. We restrict 
ourselves to the weak-coupling and singular-coupling limits, which correspond to decoherence in the energy 
eigenbasis and in the computational basis, respectively. We demonstrate that decoherence in the instantaneous 
energy eigenbasis does not necessarily detrimentally affect adiabatic quantum computation, and in particular 
that a short single-qubit T 2 time need not imply adverse consequences for the success of the quantum adiabatic 
algorithm. We further demonstrate that boundary cancellation methods, designed to improve the fidelity of adi¬ 
abatic quantum computing in the closed system setting, remain beneficial in the open system setting. To address 
the high computational cost of master equation simulations, we also demonstrate that a quantum Monte Carlo 
algorithm that explicitly accounts for a thermal bosonic bath can be used to interpolate between classical and 
quantum annealing. Our study highlights and clarifies the significantly different role played by decoherence in 
the adiabatic and circuit models of quantum computing. 

I. INTRODUCTION 

Suppose qubits with dephasing time T 2 participate in a 
computation lasting a time tf, without any error correction, 
yet T 2 ^ tf. Can this computation be quantum? It is a com¬ 
monly stated folklore position that the answer must be ‘no’. 

It is our goal in this paper to provide a theoretical basis for 
conditions under which the answer is in fact a qualified ‘yes’. 

Certainly, if one adopts the perspective of the circuit model 
of quantum computation, then the ‘no’ answer is fully justi¬ 
fied im. However, what is true in the circuit model does not 
necessarily apply directly to other models of quantum com¬ 
putation, in particular the adiabatic model ||2l, in spite of the 
fact that the two models are computationally equivalent SSI. 

Thus, a commonly held belief, that a short single-qubit de¬ 
phasing time necessarily implies quantum computational fail¬ 
ure, should not be applied without first carefully specifying 
the computational model. 

We are motivated to revisit the question of the role of deco¬ 
herence in adiabatic quantum computation (AQC) and quan¬ 
tum annealing (QA) Ib HTSll by recent experiments involving 
increasingly larger numbers of niobium superconducting flux 
qubits using programmable quantum annealing devices built 
by D-Wave llT9ll^ . For such flux qubits, the T 2 time can 
at present range from tens of nanoseconds to a few hundred 
nanoseconds EDEa, yet the computation lasts on the order 
of microseconds to milliseconds. If the qubits have all deco¬ 
hered long before the computation is over, how can this be rec¬ 
onciled with evidence that the D-Wave devices perform quan- 
7 In essence, the answer boils down to 

• The computation takes place in (or close to) the ground 
state, and 


turn annealing Il23l - l32ll 
two key points; 


• Decoherence takes place in (or close to) the instanta¬ 
neous energy eigenbasis. 

Using an adiabatic quantum master equation framework m, 
we shall explain how this leads to a very different behavior 
of adiabatic quantum computation in the presence of deco¬ 
herence than what can be expected by direct analogy from 
the circuit model. The master equation provides a consistent 
framework for the analysis of decoherence in adiabatic quan¬ 
tum computation, which will help to clear up possible mis¬ 
conceptions arising from analogies drawn too closely with the 
circuit model. 

The structure of this paper is as follows. In Sec. we 
give a brief qualitative discussion of the role of the different 
timescales (Ti, T 2 , and the total evolution time) in the cir¬ 
cuit and adiabatic quantum computing models. In Sec. Ill we 


summarize the adiabatic quantum master equation derived in 
Ref lfT4ll . to set up the main tools used in this work. In par¬ 
ticular, we distinguish between the weak and strong coupling 
limits (WCL and SCL, respectively), a distinction that gives 
rise to two very different master equations, and turns out to be 
crucial in understanding the role of decoherence in AQC. In 


Sec. IV we apply these master equations to the simple case of 


a single qubit with a time-independent Hamiltonian, and re¬ 
derive familiar results in order to establish the appearance of 
the Ti and T 2 times. We venture into new territory in Sec.fV] 
where we apply the WCL and SCL master equations to the 
case of a single-qubit coupled to a time-dependent Hamilto¬ 
nian. We consider both the adiabatic and non-adiabatic cases, 
and demonstrate explicitly that the role of decoherence is very 
different depending on whether the WCL or SCL applies. In 
particular, in the WCL thermally assisted AQC ifTSl can take 
place in the adiabatic limit, while in the SCL the population is 
distributed equally in the long time limit between the ground 
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and excited state, so that AQC becomes impossible. We also 
demonstrate that the adiabatic limit is not optimal, in the sense 
that the ground state population is maximized at a total evolu¬ 
tion time that can be much shorter than the adiabatic timescale 
and is determined by the bath spectral density. We digress 
in Sec. |Vl] to analyze the role of the interpolating function 
between the initial and final Hamiltonians on enhancing the 
probability of finding the ground state. We show that a strat¬ 
egy of imposing smooth boundary conditions developed for 
closed-system AQC in Ref. If33l (see also 041 ) has a benefi¬ 
cial effect in the open-system setting as well, though the effect 
is milder. We then come to the analysis of AQC in the multi¬ 
qubit case, in Sec. VII Here is where we address the main 
question motivating this work, namely, the role of the single¬ 
qubit T 2 on the success probability of AQC. We again demon¬ 
strate that the answer depends drastically on whether the WCL 
or SCL applies, with the former’s ground state population not 
exhibiting a dependence on the single-qubit T 2 , under reason¬ 
able non-degeneracy assumptions. In Sec. VIII we address a 
deficiency of our methodology, namely the fact that master 
equations are limited to a relatively small number of qubits 
(currently < 15), and the fact that there exists an intermediate 
coupling regime between the WCL and the SCL. We address 
this by considering a quantum Monte Carlo approach that ex¬ 
plicitly accounts for a bosonic bath, and demonstrate its util¬ 
ity in interpolating between the quantum and classical regimes 
using a quantum annealing problem that has been studied in 
the context of quantumness tests of the D-Wave device. We 
conclude in Sec. |IX| and provide additional technical details 
in the Appendix. 


II. TIMESCALES AND DECOHERENCE IN THE 
CIRCUIT VS THE ADIABATIC MODEL 


The interaction between a quantum system and its environ¬ 
ment is responsible for decoherence, which can undermine 
the efficiency of the quantum computation or even render it 
useless in the sense that it can be efficiently simulated by a 
classical computer m. Two distinct decoherence time-scales 
are usually singled out. The first is the loss of phase coher¬ 
ence between states, an elastic (energy conserving) process 
with a time-scale often referred to as T 2 . The second is the 
thermal equilibration time, which involves energy exchange 
with the thermal environment, with a time-scale often referred 
to as Ti. In the simplest case of time-independent Marko¬ 
vian dynamics, where a system can be described by a Lind- 
blad equation Il35l43^ . there is a fundamental relation relating 
these two times scales, T 2 < 2Ti ll^ , but often it is the case 
thatTa < Ti. 

In the quantum circuit paradigm, a quantum computation 
is a sequence of unitary operations (quantum gates) acting on 
one or more qubits. Quantum information is stored not only 
in the strings of O’s and I’s representing the state of the en¬ 
tire system in the computational basis, but also in the relative 
phase between superpositions of computational basis states. 
While models of circuit model quantum computation exist 
which allow for a high degree of decoherence and still enable 


a speedup over classical algorithms iolEIl, it is clear that 
in general T 2 represents an upper limit on the time it takes to 
perform a circuit model quantum computation, in the absence 
of quantum error correction 

In the adiabatic quantum computing paradigm 12, or in 
quantum annealing ll43l l44l . a computation is performed by 
evolving a system using a time-dependent Hamiltonian, with 
the final ground state encoding the solution of the computa¬ 
tional problem. Thus the energy eigenbasis replaces the com¬ 
putational basis as the relevant basis for the computation, ex¬ 
cept possibly at the end, when the system may be measured in 
the computational basis, as is typical in adiabatic quantum op¬ 
timization MM- In this setting, phase coherence between 
energy eigenstates is irrelevant. Moreover, if a near-optimal 
solution is acceptable then a computation which ends in a suf¬ 
ficiently low-lying excited state can be good enough and the 
condition that the computation terminates in the ground state 
can be relaxed ISl (though it follows from the PCP theo¬ 
rem that this does not necessarily change the complexity class 
BSl ). Therefore, as long as the final measurement can clearly 
distinguish energy eigenstates then decoherence between en¬ 
ergy states is harmless. For this reason, unlike the circuit 
model, adiabatic quantum computation is believed to exhibit 
a degree of inherent robustness to decoherence IMIol. 

Let us now give a heuristic argument why the role of the to¬ 
tal computation (i.e., evolution) time Q is quite different if one 
compares closed to open-system adiabatic quantum computa¬ 
tion. In the closed-system setting, the only relevant time-scale 
is the condition that the evolution be sufficiently adiabatic, 
i.e.. Tad ~ 1/Aini„, where Ami„ is the minimum energy gap 
between the instantaneous ground state and all excited states 
that do not become part of the ground subspace. In the open- 
system setting, if phase decoherence does occur only in the 
energy eigenbasis, the remaining relevant time-scales deter¬ 
mining the efficiency of the computation are tf, the relax¬ 
ation time Ti, and the time-scale associated with the (closed) 
system evolution being sufficiently adiabatic Tad. The in¬ 
terplay between these time-scales is non-monotonic and cer¬ 
tainly more complicated than in the closed-system setting. In 
the latter, setting the heuristic adiabatic condition, tf ^ Tad, 
guarantees—^by suppressing non-adiabatic transitions—that 
the final state reached has high overlap with the ground state 
of the final Hamiltonian lf3^ 14914521 . However, in the pres¬ 
ence of a thermal bath, even if tf ^ 1/Amin there may still 
be significant loss of population from the ground state due 
to thermal processes. The reason is that thermal excitation 
rates from the ground state typically grow as the gap shrinks 
and as Q grows, allowing the system to thermally relax into a 
Gibbs state that may have a significant population in low-lying 
energy eigenstates. Such thermal relaxation can adversely 
impact the efficiency of the adiabatic quantum computation. 
Therefore, we find that the role of Q, which was unambiguous 
in the closed-system case, becomes ambiguous when thermal 
processes are considered, and in general one can expect an 
optimal value of Q that is problem dependent fT] [l4] l5Tl455l . 
Indeed, we shall demonstrate this in Sec. [V| below for a spe¬ 
cific model involving a qubit coupled to a bosonic bath with 
an Ohmic spectral density. 
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III. MASTER EQUATIONS 

The considerations above can be made rigorous by analyz¬ 
ing a system evolving in the presence of a thermal bath that is 
described in terms of an adiabatic master equation with time- 
dependent Lindblad operators m. Such a master equation 
description has the attractive feature that it guarantees the pos¬ 
itivity of the density matrix at all times, but naturally requires 
certain assumptions and approximations. 

Consider a time-dependent system Hamiltonian 

Hs{t) \£ait)) = Sait) \£a{t)) , ( 1 ) 

where the states {|ea(i))} are the instantaneous energy eigen¬ 
states and the aforementioned gap is 


Amin = min(£a(i) - eo(i)) > 0, (2) 

a,r 

where |eo(f)) is the instantaneous ground state and |ea(f)) 
(a > 1) are the excited states. The condition Amin > 0 en¬ 
sures that only excited states that do not eventually become 
part of the ground subspace are considered. 

Next consider a generic system-bath Hamiltonian, 

Hit) = Hsit)(»tB + ts^HB + Hi (3a) 
Hi = Ba , (3b) 

Ol 

where and in the interaction Hamiltonian are, respec¬ 
tively, dimensionless Hermitian system and bath operators and 
g is the system-bath coupling strength. An adiabatic mas¬ 
ter equation in Lindblad form ll35l for the system’s evolution 
can be derived in the weak coupling limit (WCL)—where Hs 
dominates Hs b in the sense of Eq. ( |^ below—^by invoking 
the standard Born-Markov and rotating wave approximations, 
along with an adiabatic approximation ca. 

Consider the bath correlation functions (we set H = 1 
henceforth): 

Ba^^it) = . (4) 

The characteristic decay time tb is then defined via 

I (Ba^flit)) I = \TT[pBBa,0it)]\ ^ e-*/"-, (5) 

where ps is the initial state of the bath. Note that this expo¬ 
nential decay is not guaranteed but simply assumed here in 
order to extract the timescale tb- 
Now assume: 

g'^TB < Amin (weak coupling) (6a) 

gTB ^ 1 (Markov approximation) (6b) 

^ < min{A^i„,T-2}, (6c) 

tf 

where h = maXsg[o,i];a,b I (£a(s)| 9siT(s) |£b(s)) | estimates 
the rate of change of the Hamiltonian. Inequality ( [6cl i com¬ 
bines the heuristic adiabatic approximation with the condi¬ 


tion that the instantaneous energy eigenbasis should be slowly 
varying on the timescale of the bath M- 

Provided these conditions are satisfied the quantum adia¬ 
batic master equation takes the generic form HI: 

^pit) = -i [Hsit) + HLsit), pit)] + /:wcL[p(t)] 

(7a) 

t^WChlpit)] = 

OJ ^ 

-\{Li,o.it)Lfi,a,it),pit)}^ , (7b) 

where the sum over w is over the Bohr frequencies of Hs, and 
where the time-dependent Lindblad operators are 


La,Loit) = {Sait)\Aa |£&(f)) |£a (f)) (^b (f) | • 

u>=eb{t) — ea(t) 


( 8 ) 


The decay rates 

/ OO 

dt (Basit)) (9) 

-OO 


are Fourier transforms of the bath correlation function form¬ 
ing a positive matrix 7 ( 0 ;) whose elements satisfy the KMS 

conditiorQ 


7a/3(-w) = e '^“7 pq,(w) , (10) 

where (3 is the inverse temperature, and 

a/3 

is a Lamb shift term, where 

Sapioj) = [ dui'-fapiio')V ( -, ( 12 ) 

7-00 \LV-U} J 

with V denoting the Cauchy principal value. 

The Lindblad master equation in Eq. Q can be thought of 
as a natural generalization of the time-independent case 138], 
where the Hamiltonian varies sufficiently slowly relative to 
the bath so that at any instant in time we simply have a copy 
of the time-independent master equation. In fact, we can re¬ 
cover the time-independent result by replacing Hsit) by a 
time-independent system Hamiltonian. We will show in sub¬ 
sequent examples that, because of the form of the Lindblad 
operators in Eq. decoherence occurs in the instantaneous 
energy eigenbasis. 

Another standard case is the other extreme limit, where 
Hsb dominates Hs (i.e., where inequality (| 6 all is reversed). 


^ We use a slightly different convention than Ref. 03 by including the factor 
in the definitions of 7^/3 and Sa^{u}) instead of pulling it out as in 
Eqs. (46)-(50) in Ref. [T3. 
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which is often called the singular coupling limit (SCL). The 
resulting master equation takes the form: 

^p{t) = -i [Hsit) + Hi^s, p{t)] + >CscL[/o(f)] (13a) 
Cschlpit)] = ^7a/j(0) (^App{t)Al - i {A\Afi,p{t)]^ , 

(13b) 

where now i/LS = 'Thap limit, the Lind- 

blad operators are simply the bare system operators Aa, as 
is often written down in phenomenological treatments of the 
master equation. Here decoherence occurs in the basis that 
diagonalizes these operators, when such a basis exists, as we 
will show explicitly later. 


IV. DECOHERENCE OE A SINGLE QUBIT WITH A 
TIME-INDEPENDENT HAMILTONIAN 

A. Pure dephasing 

We begin our discussion with a quick review of the sim¬ 
plest case of decoherence. We consider a single qubit with a 
time-independent Hamiltonian coupled to an arbitrary bath for 
which the conditions required for the derivation of the mas¬ 
ter equation are satisfied. The total Hamiltonian is given by 
Eq. (|^, but we shall assume that 

Hs = . Hi = ga^ ® B. (14) 

For the interaction Hamiltonian in Eq. ([T4|, there is only a 
single system operator A = = |0)(0| — |1)(1|. The eigen¬ 

states are |eo(f)) = |0) and |ei(f)) = |1). Considering the 
WCL master equation 0 and {ea{t) \ A\eb{t)) oc Sa,b, there 
is only a single Lindblad operator that is non-zero: 

L.,0 = , (15) 

as given by Eq. •Si- This follows since [Hs, Hj] = 0. There¬ 
fore, the master equation for the single qubit takes the simple 
form 

^p[t) =-i[Hs,p{t)] (16) 

- 1 - 7 ( 0 ) “ 2 ) 

where we have also used the fact that H]^s oc. 1. This equation 
can be solved analytically and gives, after expanding p{t) — 
Thi jG{o 1 } (il’ ™tl taking matrix elements of Eq. Q (in 
the computational basis, which here is equivalent to the energy 
eigenbasis): 


Pm{l) — Poo(O) — 1 — Pii{t) , (17a) 

poi(f) = + iujH)poiiSl) = Pwit) , (17b) 


where 


7^(*^) _ 

-^2 — 


1 

^ ’ 


( 18 ) 


where the ‘c’ superscript denotes the computational basis (we 
shall shortly see a second T 2 associated with the energy eigen¬ 
basis). This is the familiar pure dephasing channel, where 
only the off-diagonals elements (transverse magnetization) 

(c) 

decay with a characteristic timescale T 2 . The stronger the 
coupling to the bath [recall via Eq. 0 that 7 oc p^], the 
shorter the qubit coherence time. We note that the energy 

fc) fc) 

gap Wz plays no role in the result for , and in fact, T 2 
here is entirely determined by the spectrum of the bath corre¬ 
lation function at zero frequency. In this example there is no 
thermal relaxation (the Ti time is infinite), since the popula¬ 
tion of the energy states remains fixed, as a consequence of 
[Hs,Hi]=0. 


B. Decoherence when [77s, 77/] A 0 


Let us now replace the system Hamiltonian so that 
[Hs, Hi] ^ 0. Specifically, consider 

Hs =-^utoocr"', Hi=ga^®B. (19) 

We shall study this scenario in both the weak and singular 
coupling limits. We shall see that there is a sharp contrast 
between the two, with the WCL resulting in decoherence in 
the energy eigenbasis, while the SCL results in decoherence 
in the computational basis, just as in the previous subsection, 
when Hs and 77/ were commuting. 


1. WCL 


The energy eigenstates of Hs are jeo) = |-l-) and 

|ei) = |—) with respective eigenvalues —^ojx and ^ujx, 

where |±) = ^ (|0) ± |1)). Since |±) = |=f) the non¬ 
zero Lindblad operators are: 

Lz,Ul:„ = |+)(~| ; = |—)(-|-| . ( 20 ) 

Note that we now have a non-trivial Lamb shift term: 

Hls = '5'(a;a;)|—)(—I -I- iS'(—Wa;)|-|-)(-|-| . (21) 

Writing p{t) = {+ _} Pij\'i){j\^ taking matrix el¬ 

ements of Eq. 0. we find that the master equation for the 
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density matrix components is; 


d 

dt 

d 

dt 

_d 

dt 


P — {t) = -l{^x)p — {t) +"f{-uj^)p++{t) , (22a) 

p++{'t) = l{^x)p—{t) - -i{-uj^)p++(t) , (22b) 






-i {S{uJx) - S{-UJx) + LOx) 

p-+it) . (22c) 


where we have used the KMS condition ( [TOl i to simplify the 
expressions. These equations can be solved analytically to 
give: 


p_+(t) = , (23a) 

p — {t) = PGibbs(-) + [p —(0) -PGibbs(-)] ’ , 

(23b) 

P++ (i) = 1 - P—{i) > P+-(t) = P-+ (t) , (23c) 

where 

p±iiuJa;/2 

PGibbs(±) = ^ ^ (24) 


7^(^) _ 

1 7(wxKl + e-P-^) ’ ' 


2Ti^) 


(25) 


We observe three important facts about the result in 
Eq. ( |2^ . First, the decoherence occurs in the energy eigenba- 
sis, i.e., the off-diagonal components in the energy eigenbasis 
(hence the ‘e’ superscripts on Ti and T 2 ) decay exponentially 

fe) 

to zero with a timescale determined by T '2 , and this includes 
the entire contribution of the Lamb shift. Second, the popu¬ 
lations (p++,p_) approach the Gibbs state associated with 

(e) 

the Hamiltonian Hs within a timescale determined by . 

(e) (e) 

Third, the two timescales (T^ ) have a non-trivial de¬ 

pendence on the energy gap uJx- 


2. SCL 

Let us contrast this with what happens in the SCL case, 
Eq.([Tg. In this case the evolution of the off-diagonal compo¬ 
nents is most conveniently written in the computational basis: 

d 1 

■^Poo =-^^^^x{Plo-Poi) , (26a) 

d 1 

^Pii = -* 2 ^^ (Poi - Pio) , (26b) 

d 1 

^Poi = (pii - Poo) - 27(0)por , (26c) 

^Pio = i^i^x (poo - Pii) - 27 ( 0 )pio . (26d) 


This set of equations can be solved analytically for arbitrary 
initial conditions, but for brevity, let us consider the case 
where the density matrix is initially in the ground state, i.e., 
p(0) = |-|-)(-l-|. The solution is then given by: 

1 1 _ 4. 

Poo = Pit = 2 I Pot = Pio = 2 ® ^ ■ (22) 

In this case, the off-diagonal elements in the computational 

(c) 

basis decay exponentially with a timescale determined by T 2 
[Eq. ([T8|], so we have decoherence in the computational basis 
regardless of the fact that the system Hamiltonian does not 
commute with Hj. 

These simple time-independent examples anticipate what 
we shall see for their time-dependent counterparts. In the 
WCL, we expect decoherence in the energy eigenbasis, a fea¬ 
ture that does not preclude the success of an adiabatic quan¬ 
tum computation, whereas in the SCL, we expect decoherence 
in the computational basis, rendering adiabatic quantum com¬ 
putation impossible. 


V. DECOHERENCE OF A SINGLE QUBIT WITH A 
TIME-DEPENDENT HAMILTONIAN 


We consider the following time-dependent single qubit 
Hamiltonian, with a linear interpolation schedule; 

Hsit) = -^i^x (1 -s)(T^ - , (28) 

where t G [0, f/] and s = t/tj is the dimensionless time (the 
role of the interpolation schedule is studied in Sec. VI 1 . This 


Hamiltonian interpolates between the two system Hamiltoni¬ 
ans considered in Eqs. ( [T4l l and ( [T^ in the time-independent 
case above. The instantaneous energy eigenvalues are: 

s±(a) = ±2^(1^^ + ^ ^ , (29) 


and hence the instantaneous energy gap is A (5). Its minimum 
is 


Amin — 






(30) 


which is reached when s = Smin^ where 


Smin — ^ — ^zl^x ■ (31) 

1 -f i ^ 

It is convenient to define the dimensionless instantaneous gap 
\{s) = ^ = x/(l - s)2 + s2r2 , (32) 
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so that Amin = Amin/wa: = F/V1 + The Corresponding 

energy eigenstates are 


\e+{s)) = 


1 


c+(s) 


sT — X{s) 
I — s 


| 0 ) + | 1 ) 


(33a) 


|e_(s)) = —- [(sT + A(s)) |0) + (1 - s) |1)] , (33b) 

c_(s) 

with c±(s) being the appropriate normalization factors. Note 
that when s = 1 the ground state is |0). 

To identify the adiabatic limit for this case, let us use the 
heuristic adiabatic condition [Eq. 


equations of motion: 


r 


2A2(s) 


(P-+ + P+-) + [-^+( 5 )^++ - ^-{s)p—] 


(38a) 


^P+-is) = ^p*_+is) 

as as 


2X^{s) 


{p++ - P —) - + E(s)]p+_ , (38b) 


where 


max 

0<s<l 


I {e+\dsH\e-) \ 
f/A(s)2 


^/A^n 


< 1 , 


(34) 


where we have used the fact that the numerator (which equals 
rwa;/[ 2 A(s)] before maximization) and denominator are both 
respectively maximized and minimized at s = Smin- Rewrit¬ 
ing, this yields for the adiabatic condition 

r 

tfOJx 3> — 3 — or, equivalently (35a) 

^ (r + , (35b) 


where the second form of the inequality emphasizes the sym¬ 
metry between and . 


-7^±(s) = f/C^(s) 7 (±A(s)) (39a) 

n{s) = tf [A(s) + {5(A(s)) - ^(-A(s))} C^(s)] 

(39b) 


E(s) = tf 



\ [7(A(s)) -f 7 (-A(s))] C^(s) 


(39c) 


The term of Eq. ( |38b| ) involving S(s) gives rise to the ex¬ 
ponential decay of the off-diagonal elements of the density 
matrix in the instantaneous energy eigenbasis. The corre¬ 
sponding decoherence time [f//S(s) in Eq. ( |39cl l] interpo¬ 
lates between the and times given in Eqs. ( [T8| ) and 
( |^ , respectively: fy/E(0) = 2 /[ 7 (u; 2 ,)(l -f 
[where we used the KMS condition ( [TOl l] and fy/E(l) = 

1 /[ 27 ( 0 )] = T 2 ‘^\ This is to be expected since as mentioned 
above, Hs{t) [Eq. ( [28] l] interpolates between the correspond¬ 
ing two system Hamiltonians. 


A. WCL 


In the WCL, the Lindblad operators take the form: 


1. Solution in the adiabatic limit 


= yTA (I^-('S))(£-(s)I - k-r(s))(^-r(s)l) (36a) 
A(s) 

L±\is) = Cis)\e^{s)){e±{s)\, (36b) 

where 

C(.) = . (37) 

(We have dropped the index a from the Lindblad opera¬ 
tors since there is only a single qubit.) Note that 0 < 
C(s) < 1 and that these Lindblad operators interpolate 

between the operators in Eqs. (ED and ( |20l l. Denoting 
Pif{s) = (ej(s)lp(s) |£^(s)) (with e {+,-}), and us¬ 
ing (e+l ds |e-) = — (e- I ds |e+) = —T/2A(s)2, we have as 


The first summands in Eq. ( [38] l (proportional to E/[2A^ (s)]) 
are purely due to the evolving instantaneous energy, and are a 
factor of f/ smaller than the remaining terms. Therefore, for 
sufficiently large tj (i.e., when the adiabatic condition ( |35a| ) 
is well satisfied), these terms can be neglected. It is simpler to 
analytically solve the dynamical equations ( |38l l in this limit, 
since it decouples the diagonal and off-diagonal elements. Us¬ 
ing the KMS condition ( [T0| ) again to relate J^+(s) and jF-{s) 
we can rewrite these equations in the adiabatic limit as: 


^p__(s) = T+{s) [1 - (1 + p-{s) 


, (40a) 


ds 


p+-{s) = - [in{s) + E(s)] p+-(s) . 


(40b) 
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These equations have solutions given by: 



P+- (s) = exp 


ds' + S(s')) 


P++(s) = 1 - P—{s) , 
P- + is) = P+-(s) • 


P+-(0) , 

(41b) 

(41c) 

(41d) 


As is already clear from Eq. ( |40a| ), any deviation in the popu¬ 
lation of the instantaneous ground state, i.e., p _, in the adi¬ 

abatic limit is purely due to a non-zero J 74 .(s), which in turn 
requires a non-zero 7 (A(s)), i.e., a “resonant” thermal excita¬ 
tion. This means that the rate of population loss is expected to 
be non-monotonic in the instantaneous energy gap A(s), i.e., 
when the goal is to prevent population loss from the ground 
state, increasing the gap can do damage before it starts to 
help. The details of this depend on the decay rate 7 (w) via the 
bath correlation function B{f), as is clear from Eq. (see 
also Eig.[^. Also noteworthy is that the Lamb shift S{A{s)) 
affects only the off-diagonal elements, through r 2 (s). 

Equation ( |40b| l highlights that, even in the time-dependent 
case, the decoherence occurs in the instantaneous energy 
eigenbasis. In the adiabatic limit in which Eq. ( |40| is derived, 
the dynamics of the phase coherence between energy eigen¬ 
states completely decouples from that of the energy state pop¬ 
ulations, such that its exponential decay [at a rate determined 
by E(s), Eq. ( |39c| l] does not affect the evolution of the energy 

state populations. If p-|_is initially zero, e.g., if the system is 

initialized purely in the ground state, then no coherence with 
the excited state is ever generated. Therefore, we find that 
in the WCL, even in the presence of fast dephasing, an adia¬ 
batic computation is possible. Similarly, Eq. ( |40a| i highlights 
that the thermal excitation/relaxation processes occur in the 
instantaneous energy eigenbasis since the dynamics is entirely 
determined by the population in the energy eigenstates 


2. Solution without the adiabatic limit 


Let us now return to the full dynamical equations ( [38] l with¬ 
out taking the adiabatic limit. Since T, and E depend on 
tf, there is a non-trivial dependence on tf for the final ground 
state population. To explore this dependence we solve the dy¬ 
namical equations ( [38l l numerically. We assume the bath is in 
a thermal state with an Ohmic spectral density, i.e., 

7 ( 0 ;) = 2TTr,g , (42) 


y M 



FIG. 1. Ohmic spectral density [Eq. @] with 1/P = 

2.23GHz, tJc = SttGHz. 


where Wc is a high-frequency cutoff and 7 is a positive con¬ 
stant with dimensions of time squared arising in the specifi¬ 
cation of the Ohmic spectral function. This spectral density 
is depicted in Fig. for typical parameters used in our nu¬ 
merical calculations. When the cutoff is the largest energy 
scale, specifically when Wc ^ 1/P, the bath correlation time 
[Eq. (|^] can be shown to be tb = p/(27r) |[l4l . Thus the 
validity condition (| 6 c|l becomes, in addition to Eq. (|35a|l: 


f/» 



or, equivalently 


tf^/uJfuTz » - 



(43a) 

(43b) 


Examples of evolutions satisfying Eq. (gg are shown in 
Fig. 0 for both the closed system case and the open sys¬ 
tem case for increasing values of tf. It shows that oscilla¬ 
tions of the instantaneous ground state population due to non- 
adiabatic unitary dynamics are damped out as the total evo¬ 
lution increases, and that some ground state population is re¬ 
covered due to the thermal relaxation for large enough tf. The 
dependence of the final population on tf is shown in Fig. 
where we explicitly see the effect of the various timescales in 
our problem. For very short evolution times (where the adi¬ 
abatic condition is not satisfied, i.e., 2tfUJxX^^^/T < 1 , the 
evolution is highly non-adiabatic, and the final ground state 
probability is close to 1/2 (see the left inset of Fig. |^. How¬ 
ever, note that in this regime we may not entirely trust the 
ME to be a reliable approximation for the dynamics since the 
condition h/tf < [Eq. @] requiring that the system 
evolves much more slowly than the time-scale of the bath is 
not necessarily satisfied. Furthermore, since the evolution is 
so short that 7^+(s) < tfj(A(0)) = tfj(ccx) <C 1 [recall 
Eq. P 8 a[ )], thermal effects are small because they have insuf¬ 
ficient time to act. 

As we increase tf, the evolution becomes more and more 
adiabatic and the ground state probability peaks close to 1 
around 2tfUJx\%,{xL/^ — 1’ '^^e adiabatic condition P5a| i (see 
first inset of Fig. |^. However, as we continue to increase tf, 
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P-is) 




p-{s) 



P—is) 



FIG. 2. Time-dependence of the ground state population in the WCL for a quhit evolving adiahatically subject to the Hamiltonian ( |28t with 
p(0) = |£-(0)) (e_(0)|. In (a) we depict the case of closed system evolution with tfOJx = 10%/2, where the oscillation is entirely due to 
unitary non-adiabatic transitions. In (b)-(d) the system (with tfUJx = lOv^. tfOOx = 5 x 10^, tfuix = 5 x 10^ respectively) is coupled 
to an Ohmic bath with 'y(uj) as in Eq. l |42^ , and the instantaneous ground state population becomes gradually more damped as tf increases 
[the insets in (c) and (d) zoom in on the vertical axis]. The instantaneous ground state population initially decreases rapidly due to thermal 
excitation and then recovers somewhat due to thermal relaxation. In (b) the oscillation is a damped version of what is seen in (a). In (c) the 
oscillation is completely damped by the bath, and a small population recovery is seen towards the end of the evolution. This recovery grows 
at even larger tf as seen in (d), indicating that this is due to thermal relaxation (this is shown in more detail in Fig.|^. Parameters are chosen 
to satisfy the various inequalities: ujx = Uz = 1 GHz so T = 1 and Amin = 1 / V2 and Eq. j35at is satisfied sincer/2Amin = \/i^tfUJx- 
Eor the Ohmic spectral density we chose = 10“"^, 1/|3 = 2.23GHz < uj^ = SttGHz (as in Fig.[^ so that Eq. j43[ ) applies and is satisfied 

since ^ ^ < i < 


rather than observing that the system remains in its ground 
state, thermal excitations increase; this removes significant 
population from the ground state, which actually drops below 
its thermal equilibrium probability distribution. An example 
of this is shown in Fig. 2(c) where the ground state population 
decreases over the first half of the evolution. As we continue 
to increase tf, thermal relaxation allows the system to relax 
to its thermal equilibrium probability distribution. For an ex¬ 
ample of the effect of thermal relaxation, see the increase in 
ground state probability during the second half of the evolu¬ 
tion in Fig. |2(d)| The increase in the population as we increase 
tf is shown in the second inset of Fig. 


Thus thermal excitations can have a significant detrimental 
impact on the ground state population, and hence the success 
probability of an adiabatic quantum computation. As can be 
seen in Fig. there is an optimum value for t / that maxi¬ 


mizes the ground state probability. This value balances the 
adiabaticity of the evolution against the time allowed for ther¬ 
mal processes to occur. We illustrate the dependence of this 
optimal time on the bath strength in Fig. where we see that 
as the system-bath strength increases in value, such that ther¬ 
mal processes occur more rapidly, the optimal evolution time 
decreases. This shows that in open system adiabatic quan¬ 
tum computation it can be advantageous to stop the compu¬ 
tation early, before the ground state probability starts to dip 
due to thermal excitations, in agreement with earlier findings 

EiHisiia. 
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FIG. 3. Final WCL ground state population (Pqs = p _ i^f)) a 

function of evolution time [in dimensionless units corresponding to 
the adiabatic condition j35aH for the model specified in Fig. with 
the same initial condition. The left inset zooms in on the short total¬ 
time evolution, the right inset shows the long total-time evolution. 
The maximum in the ground state probability (left inset) is the opti¬ 
mal evolution time and is seen to occur at 2tfL0j:X^i^/r « 7, i.e., 
increasing tf much above the heuristic adiabatic condition does not 
help. The reason is that this maximum is a balance between maximiz¬ 
ing adiabaticity while minimizing thermal excitations. The right in¬ 
set shows that for very long evolution times thermal relaxation repop¬ 
ulates the ground state, which eventually settles on its Gibbs distri¬ 
bution value. The condition j43| l [expressed as (^) — 1] 

is satisfied already at 2tfUJxXmin/^ ~ 10“®. 



(DcP 

FIG. 4. Behavior of the the “adiabatic parameter” defined in 
Eq. j35a[ l at the optimal evolution time fj' (that maximizes the 
ground state probability) on the system-bath coupling. As the cou¬ 
pling strength increases, the optimal evolution time becomes smaller, 
and the evolution becomes less adiabatic. Simulation parameters: 

= IGHz, r = 1, 1/|3 = 2.23GHz, cuc = SttGHz. Note 
that the adiabatic parameter never dips below 1 (which would vi¬ 
olate the heuristic adiabatic condition): fj' no longer exists past 
the last data point shown because beyond that point the highest Pqs 
achievable is the thermal one, which cannot be exceeded no matter 
how large tf becomes. The smallest system-bath coupling plotted is 
^ = 5.67 X 10“®. 

UlcP 


B. SCL 


We can perform a similar analysis in the case of the SCL, 
Eq. ([^. We find that the dynamical equations are given by: 

d d * \ / \ A ■s 

^Poo = ~ — s) (poi ~ Pio) j (44a) 

d _ _ tf 

^P01-^Pl0--^P01 + 

-^2 

[(1 - s) (pii - Poo) + 2srpoi]. (44b) 


We plot their numerical solution in Figs.|^and|^ With the ini¬ 
tial condition being a fully populated ground state, we observe 
behavior that depends strongly on the various timescales. For 
tf < = [ 27 ( 0 )]“^ but tf » r/(a; 2 ;A^i„) [the adia¬ 

batic condition P5a[ l], the decoherence in the computational 
basis does not have enough time to disrupt the adiabatic evo¬ 
lution, and the system can perform an adiabatic computation. 
This is illustrated in Fig. |5(a)| where the population in the 
ground state at t he en d of the evolution is almost 1. But when 

tf > 4"^ [Fig. 


5(b) I, even if the adiabatic condition is well 


satisfied, the off-diagonal entries of the density matrix start to 
decay exponentially and we begin to see a significant loss in 
the ground state population at the end of the evolution. For 
tf T 2 , the off-diagonal terms rapidly decay to zero, and 
the evolution results in a state that resembles the maximally 
mixed state [Fig. 5(c) |. In this case the evolution is very well 


approximated by the results of the time-independent case in 
Eq. ( |Z7] l. The dependence on the final ground state probability 
is shown in Fig.|^ and is seen to rapidly tend to the maximally 
mixed state value of 1 / 2 . 


C. WCL vs SCL 


Concluding this study of the single qubit case, we observe 
some important differences between the WCL and the SCL. 
Most notably, as is apparent from comparing Figs. [^and|^ 
while superficially the overall behavior of the ground state 
population appears qualitatively similar (especially for short 
tf), the key difference is that in the WCL the final ground 
state population settles on a finite-temperature thermal equi¬ 
librium value, while in the SCL it approaches the infinite tem¬ 
perature (maximally mixed state) thermal equilibrium value of 
1/2. We shall see these conclusions reinforced below, when 
we discuss the multi-qubit case. The implication is that adi¬ 
abatic quantum computation is possible in the WCL (since 
there is a non-vanishing probability for the system to end up 
in its ground state), while it is hopeless in the SCL where 
the final equilibrium value is 1/2^, where N is the system 
size (number of qubits). The explanation for these conclu¬ 
sions is that in the case of the WCL the relevant timescales 
are the adiabatic and thermal relaxation timescales. The de¬ 
phasing timescale does not matter since dephasing in the en¬ 
ergy eigenbasis is not detrimental to the process of finding the 
final ground state. In contrast, in the case of the SCL the rel¬ 
evant timescales are the adiabatic and dephasing timescales. 
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Poo(-s) 


Poo(s) 



Poo(s) 



FIG. 5. Time-dependence of the ground state population in the SCL for a quhit evolving adiabatically subject to the Hamiltonian l |28| l and 
coupled to an Ohmic bath with 7 (tp) as in Eq. yielding 7 ( 0 ) = 21X7)0^/fi. The initial state is p(0) = |e_(0)) (£_(0)| = |+){+|- 
Depicted is the numerical solution of Eq. j44| l. In (a) uj^tf ~ 10^ such that T tf so that the evolution is adiabatic 

and unaffected by decoherence, and the ground state population increases to 1 in the computational basis. In (b) uj^tf = 10® such that 
r/(tPa;Amin) tf < resulting in a significant loss of ground state population. In (c) uixtf = 10"^ such that tf ^ and the ground 
state population rapidly settles on 1/2. The temperature and spectral density parameters are as in Fig.|^ 



'2u}xtf AiYiin / r 


FIG. 6 . Final SCL ground state population [Pgs = poo(f/)] as a 
function of evolution time. The left inset zooms in on the short 
total-time evolution, the right inset shows the long total-time evo¬ 
lution. For these long total-time evolutions, the system state rapidly 
becomes maximally mixed. The parameters are as in Fig.[^ 


system results have also been reported ll59l . In particular, 
Ref. 13^ showed that in the closed system case, the deviation 
from the final ground state can be made arbitrarily small at a 
fixed tf (large enough to satisfy the adiabatic condition with 
respect to the minimum energy gap) by choosing a family of 
analytic interpolation functions with an increasing number k 
of vanishing derivatives at the boundaries f = 0 and t = tf. 
In a nutshell, each such vanishing derivative cancels another 
boundary term in an integration-by-parts version of adiabatic 
perturbation theory. In Ref. ||58]| an explicit example of such 
an interpolation function was provided, again for the closed 
system setting, which we revisit here in the open system set¬ 
ting. Specifically, Ref. Il58ll proposed the interpolating func¬ 
tion to be given by the regularized incomplete beta function 

Sk{s): 


Bs{l -F /c, 1 + fc) 

Bi (1 "F fc, 1 -F fc) 


(45) 


Since the dephasing is occurring in the computational basis, 
this destroys the coherence of the energy eigenstates. 


VI. BOUNDARY CANCELLATION METHOD 

So far, we have only considered a linear interpolating func¬ 
tion between the initial and final Hamiltonian. However, it 
is well known that optimization of the interpolation function 
can have important consequences. For example, it is only af¬ 
ter such an optimization, based on a time-local adiabatic con¬ 
dition, that the quadratic speedup of Grover’s algorithm was 
shown to be realizable in AQC ll56l . Subsequently various 
studies have shown that the interpolating function can be opti¬ 
mized to improve the performance of the adiabatic algorithm 
in the closed system setting l34l l57l l58ll . Related open 


Okis) 



FIG. 7. The regularized incomplete beta functions [Eq. 65| >] for k = 
0 (blue-solid, linear), 1 (green-dashed), 2 (red-dotted), 5 (orange- 
dot-dashed), 10 (purple-long dashed, steepest rise). 
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log(l-PGs) 



log(l-PGs) 


I0g(l-/>GS) 




Ulftf 


FIG. 8. Log base 10 of the ground state error (1 minus the ground state probability) as a function of the total evolution time tf (in units of uJx) 
for different regularized incomplete beta function [Eq. l |45[ >] values of A: = 0 (blue, solid), fc = 1 (green, dashed), k = 2 (red, dotted). The 
three cases shown correspond to (a) the closed system case, 'qg^ = 0, (b) weak system-bath coupling, qg^ = 10“®, (c) somewhat stronger 
system-bath coupling, qg^ = 10~'*. The efficacy of the boundary cancellation method decreases as the system-bath coupling strength grows. 


where Bs{a, b) = dy y°' ^(1 — y)^ ^ with 3?(a), ^{b) > 
0 and |s| < 1- A sample of the behavior of the family 
parametrized by k (the number of vanishing derivatives at the 
boundaries) is shown in Fig. We again consider a single 
qubit in the WCL with the system Hamiltonian 

Hs{t) = - [1 - 0{s)] - ^6>(s)a;zcr^ , (46) 

coupled to an Ohmic bath, as in Sec. |V A| Figure displays 
the behavior of the ground state probability as we increase k 
starting from zero (linear interpolation). We observe that even 
in the open system case increasing k leads to an improvement 
in the ground state probability in the adiabatic regime, though 
the improvement rapidly saturates as k grows. Another attrac¬ 
tive feature of increasing k, visible in Fig. is that it sup¬ 
presses the oscillations in the ground state probability. 

It is interesting to check whether the improvement in the 
ground state probability is related to a change in the spectrum 
induced by the choice of interpolation function. As shown in 
Fig. 1^ the closed-system minimum energy gap remains un¬ 
changed by changing fc, but increasing fc results in the gap 


A/oJx 



FIG. 9. The energy gap between the ground state and first excited 
state for the annealing schedules parametrized by fc = 0,1, 2, 5,10 
(blue-solid, green-dashed, red-dotted, orange-dot-dashed, purple- 
long dashed). The minimum gap remains unchanged, but the gap 
remains large for longer as fc increases. 


remaining large for a longer time. This, in turn, is corre¬ 
lated with a reduction in the amount of thermal excitations, 


as shown in Fig. 10 


Thus, the boundary cancellation method has a positive, al¬ 
beit mild effect, even for open systems. We next turn our at¬ 
tention to the multi-qubit case. 


VII. THE MULTI-QUBIT CASE 

We now extend our discussion to the multi-qubit case, 
where the differences between the SCL and the WCL become 
most prominent. In particular, we shall address the key ques¬ 
tion posed in the introduction: how exactly does the single¬ 
qubit T 2 enter, and how detrimental is a short T 2 for successful 
AQC? We shall see that the answer depends dramatically on 
whether the SCL or WCL applies. Namely, we shall demon¬ 
strate that in the SCL with independent baths, the decoher¬ 
ence time can scale with the system size and hence rapidly 
destroy the coherence of energy eigenstates. This renders any 
computation in AQC effectively impossible in the SCL as the 
system size grows, without extensive error correction. On the 
other hand, we shall demonstrate that even in the multi-qubit 
case, AQC remains possible in the WCL limit, with the dom¬ 
inant source of error being thermal relaxation to the finite- 
temperature thermal state, but with the single-qubit T 2 not 


Iog(l-PGS(0) t 



FIG. 10. Evolution of the instantaneous ground state error for dif¬ 
ferent regularized incomplete beta function values of fc = 0 (blue- 
solid), fc = 1 (green-dashed), fc = 2 (red-dotted). We set qg^ = 
10“^, and = 100. 
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playing an important role. 


A. SCL 


Comparing Figs. and it might appear that the results 
at the optimal evolution time are the same for the WCL and 
SCL case. This is a result we expect to hold only for the sin¬ 
gle qubit case. As we increase the number of qubits N, the 
decoherence time scales with N in such a way that we would 
not expect AQC to be possible at all in the case of the SCL. To 
see this, we start from the general SCL master equation ( [T3| l 
and assume from now on that the the system operators {Aa} 
are single-qubit operators, i.e., that the index a € {1,..., N} 
enumerates the qubits. Our pessimistic conclusions about the 
SCL will not be improved by considering the case of gen¬ 
eral decoherence (e.g., with different Pauli operators acting 
on a given qubit), so this assumption will serve to illustrate 
the limitations of AQC in the SCL. Moreover, we may as¬ 
sume without loss of generality that the {Aa} are Hermitian. 
In this case the {Aa} operators commute and there exists a 
mutual diagonalizing basis {|a)}. In this basis we can write 
Aa\a) = Aaa |a), where Aaa = {a\Aa\a) (the ath eigen¬ 
value of Aa), and pat = (a| pit) \b). We consider the contri¬ 
bution of the dissipative part and obtain 


{a\£scL[p{t)] \b) 

— ^ ^ '~fo£p (Q) '^-^cx.a-^Pa 

— ^ ^ To;/3(0) (^aa 
a.,/3 


(47a) 

'^AabApj^ Pah 

Apb) Pab (47b) 


where we have used 7 a/ 3 ( 0 ) = 7 / 3 q( 0 ) [which follows from 
the KMS condition ([TOll]. Let us now consider two opposite 
extremes of decoherence. 


I. Independent Decoherence 

In the case of identical, independent baths, we have 

7 a/3(0) = (5a/37(0)- Thus 

(a| Cschlpit)] \b) = -TabPab , (48) 

where 

1 ^ 

[(4laa - Aab)/2f (49) 

is the decay rate of pab in the diagonalizing basis of Aa, and 
where = l/[ 27 ( 0 )] is the single-qubit dephasing time 
[Eq. ^]. 

When the {Aa} are Pauli operators the eigenvalues are ±1, 
so that Aaa = (—1)““, where Oq S {0,1}. Thus \Aaa — 
Aab\ = ‘2,da,ab where da,ab € (0,1} is the Hamming distance 


between | a) and | b) on the ath qubit, and hence 

= (50) 

-^2 OL ^2 

Combining Eqs. ( |49| ) and ( [50l l we see that Vaa = 0, i.e., the 
dissipative part does not affect the populations directly, and all 
off-diagonal elements decay exponentially with a rate Tab > 
0, which can be as high as N times the single-qubit dephasing 
rate. 

While we have not explicitly demonstrated that the popu¬ 
lations equalize in this setting (as we saw in the single-qubit 
case in Sec. |VB| i, it is still clear that no useful AQC can take 
place: the ground state of iTs(f) will, in general, be a coherent 
superposition of the complete set of eigenstates of the {Aa} 
operators, and we have demonstrated that this superposition 
decays on multiple timescales, varying from the single-qubit 
dephasing time to N times this timescale. Thus, to be 
able to perform useful AQC in the SCL under the indepen¬ 
dent decoherence model, one must invoke some form of error 
correction, suppression, or avoidance ll42l I60l - l6^ . 


2. Collective Decoherence 

In the case of collective decoherence there is only one sys¬ 
tem operator: A = '^a=i^a, and a single rate 7 a/ 3 ( 0 ) = 
7 ( 0 ). In the diagonalizing basis we have 

(a| £scL[p(i)] 1^) =-(7y [i^a - Ab)/2f pab ■ (51) 

The eigenvalues are — 2ha where ha S (0,1,... A^j 

is the Hamming weight of |o), with multiplicity = i^)- 
The “singlet” case Aa = 0 (which arises only when N is 
even) is the well-known decoherence-free subspace (DES) 
for collective dephasing, for which (a| £scL[p(f)] \b) = 0 , 
and we see that the DES is spanned by the (^^ 2 ) compu¬ 
tational basis states having an equal number of O’s and I’s 
B2]| . Other subspaces, defined by the degeneracy condition 
Aa = Ab, are also decoherence-free, but the latter is the 
largest such subspace. States belonging to different such sub¬ 
spaces (i.e., for which ha ^ hb) have a positive dephasing 
rate {ha — hb)^ jT^'^■ However, provided the adiabatic quan¬ 
tum computation is initialized inside a given DES, it will pro¬ 
ceed in a completely unitary manner, and its success proba¬ 
bility will be determined purely by the adiabatic condition for 
closed systems. Thus, under collective dephasing conditions 
it is possible to support adiabatic quantum computation even 
subject to the SCL. 

It is not difficult to extend this analysis to the non- 
commutative case of collective decoherence with different 
Pauli operators ll^lTOll . but collectiveness is, of course, a very 
strong condition (though it can be achieved using dynamical 
decoupling Il2l]|72l), so we shall not pursue this further here. 
As we shall see next, in the WCL case the prospects for AQC 
are significantly more favorable than in the SCL, precisely in 
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the opposite limit of an absence of any degeneracy-inducing 
symmetries. 


B. WCL 


1 . Coherence 


We showed that in the single qubit case, the WCL implies 
dephasing in the instantaneous energy eigenbasis. We now 
wish to check whether this remains true in the multi-qubit 
case, and whether AQC remains viable in the WCL. We fo¬ 
cus on the dissipative part of the WCL master equation 0 
and suppress the explicit time-dependence for notational sim¬ 
plicity: 


•^wcl[p] = 

oj ct,/3 



We consider the off-diagonal elements of this operator: 


(a| CwCl[p] I^) — 'y ^ ^ ' 'y ^ laisi^) Lf}^^^acPcdL^ci,uj,db 

u; a,/3 c,d 


-- (n 
2 


cd-^^,^,db + L^ct,LJ,acLl3,uj,cdPdb 


( 53 ) 


where now the complete set {|a)} is the energy eigenba¬ 
sis, i.e., the instantaneous eigenstates of Hs{t), and where 
La,uj,ac = (o| |c). We would like to extract the 

time from this expression. This is not possible in general, but 
as we show in detail in Appendix [A] it is possible under the 
assumption that there are no accidental symmetries, i.e., that 
the spectrum is non-degenerate, and that moreover the energy 
gaps are also non-degenerate, i.e., Se^,-e^, 6 t,-st = Sa,a'Sb,b'- 
This allows us to explicitly write the T2 time associated with 
any pair of energy eigenstates |a) and \b) (with a 7^ b): 




^ ^ 'Ta/j(Q) (-^a,aa ^a,bb) (-^/3,aa -^0,bb) 


OL,P 

1 


EE ^b'')-^a,ab' 


fdb'a 


a,/3 \b^^a 


T 'y ^ ’yajdi^^b ^a'^-^a,ba' -^f 3 ,a'b | ■ ( 54 ) 

a' ^b 


The single qubit result in Eq. ( | 39 c| i is a special case of 
Eq. (| 54 | ), as can be seen by taking a = , 5 = e+ [from 

Eq. ( |33| i] and noting that A++ = —A _= and A _^ = 

Xl-i— ■ 

In the case of identical, independent baths, we can use 


7a/3(w) = 7 (w)5c(/ 3 to further simplify this expression to: 

jlE “ if ? 

cx. \b'^a 


+ E “ Sa')\Aa, 

a'^b 


ba' 


( 55 ) 


(c) 

where we have introduced the single-qubit T2 time via 
7(0) = 1/[2T2‘^^]. Since each term in the sums is now mani¬ 
festly positive, this shows explicitly how a small single qubit 

(c) (e) 

T2 time enforces a small T2 dephasing time, but in the in¬ 
stantaneous energy eigenbasis. 


2 . Ground State Population 


(c) 

While a small T2 time causes rapid decoherence between 
energy eigenstates, it does not necessarily translate to a small 
thermalization time. To see this, starting again from Eq. ( | 5 ^ , 
we can determine the rate equations for the populations in the 
energy eigenbasis. In the absence of degeneracies (see Ap- 
pendixj^for details), we find that poo = (^ol -CwclM |eo) = 
~t‘oPoo+X]h>o’ where the sum over 6 is a relaxation term that 
repopulates the ground state, and the rate rg of depopulation 
of the ground state due to the dissipative dynamics (again for 
identical independent baths in the absence of degeneracies) is 
given by: 


ro = 7 (£a - £o)e ^ l^la.OaP ■ ( 56 ) 

a>0 a 

As is evident, the 7(0) term is absent in this sum, indicat- 
ing that the single qubit T2 does not play a detrimental role 
in depopulating the ground state. This conclusion is robust 
even in the presence of degeneracies (which we have ignored), 
since the corresponding 7(0) terms would arise due to popu¬ 
lation transfer between degenerate ground states. A problem 
would arise in that case only if a degenerate ground state be¬ 
came an excited state later in the evolution. 

We may thus conclude that AQC in the WCL is largely un- 

fc) 

affected by a small single-qubit T2 : decoherence is between 
energy eigenstates, which is harmless, and ground state de- 
population does not depend on T2 , as long the energy gap 
does not close (a scenario that is detrimental to AQC even in 
the closed system case). Ground state depopulation is pro¬ 
tected by the gap via the Boltzmann factors 
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VIII. SQA-EB: SIMULATED QUANTUM ANNEALING 
WITH AN EXPLICIT BATH 


Our master equation analysis has allowed us to study two 
limits: the WCL, where the system-bath interaction is weak 
relative to the system Hamiltonian, and the SCL, where the 
system-bath interaction is strong relative to the system Hamil¬ 
tonian. In either case, the master equation approach becomes 
computationally prohibitive when the system size becomes 
large, since in principle it scales with the square of the di¬ 
mension of the system Hilbert space. In practice we are re¬ 
stricted to simulating up to about 15 qubits in this manner. 
Furthermore, in order to interpolate between the two limits, 
an explicit treatment of the bath degrees of freedom is neces¬ 
sary. In the case of bosonic baths with dephasing interactions, 
this can be achieved by integrating out the bath degrees of 
freedom. For the case of a single qubit, analytic expression 
for the dynamics of the two-level system can then be found 
using the “non-interacting blip approximation’ ’ EH CD, and 
this method can be successful in capturing the dynamics of 
multi-qubit open system adiabatic quantum computing if the 
dynamics is effectively restricted to only two levels ED. 


-f;--f 


' -f--t- : 

FIG. 11. The 8 qubit “quantum signature” Hamiltonian F/iaing stud¬ 
ied in Ref. ED. The spins are depicted by colored disks. All 
spin-spin couplings (black lines connecting spins) are ferromag¬ 
netic with magnitude 1 and the signs of the local fields (of magni¬ 
tude 1) are indicated within the disks. We use the sign convention 

iflsing = - Ei + Ei<j 


slice, and 

J±{t) = ln(tanh(/3A(f)/A^T-)) > 0 (59) 


In this section we present a different approach, that is nu¬ 
merically efficient in the sense that its computational cost 
scales in the same manner as classical Monte Carlo methods. 
Moreover, this will allow us to probe the intermediate regime, 
between the WCL and the SCL. The price to be paid is that 
instead of simulating the dynamics we will be sampling from 
the instantaneous Gibbs distribution of the system. 

To be explicit we now restrict our attention to time- 
dependent system Hamiltonians of the form of the transverse 
Ising model. 


i 


i i<j 


(57) 

In this setting adiabatic quantum computing is also known as 
quantum annealing ||4^|751 . In simulated quantum annealing 
Il25ll76l[77l (SQA), Monte Carlo dynamics are used to sample 
from the instantaneous Gibbs state along the annealing evolu¬ 
tion. For example, in the case of discrete-time Monte Carlo, 
for each fixed time t in Eq. ( [57] ), Monte Carlo sampling of the 
thermal state associated with the transverse Ising Hamiltonian 
is done by sampling the dual classical spin system with: 




Ji 






Aj,T 


i<j 


-jmY. A^2,r/^2,T+l j 


(58) 


where j3 is the inverse temperature of the Monte Carlo simula¬ 
tion, Nr is the number of Trotter slices used along the Trotter 
direction (also referred to as the imaginary-time or time-like 
direction), yui t- denotes the ith classical spin on the rth Trotter 


is the nearest-neighbor coupling strength along the Trotter di¬ 
rection. Although this does not capture the unitary dynam¬ 
ics of the quantum system, the sampling of the instantaneous 
Gibbs state mimics the thermalization process towards the 
Gibbs state in the WCL master equation, with the advantage 
that the simulation remains efficient, so that the system size 
can be made large. However, even as the temperature is in¬ 
creased, neither the WCL master equation nor the standard 
SQA methods capture the decoherence of the qubits into clas¬ 
sical bits, unlike the SCL master equation. 

In this section we show that by incorporating an explicit 
bath in SQA, we can interpolate between SQA and classical 
simulated annealing (SA), where the qubits have fully deco¬ 
hered into classical bits that undergo bit flip updates. We call 
this method SQA-EB (EB for “explicit bath”). In addition to 
Eq. we include independent yet identical bosonic baths 
for each qubit with a dephasing system-bath interaction: 

N 

Hsb (bi^k + • (60) 

2—1 k 

Assuming the baths to have an Ohmic spectral function, we 
can analytically integrate out the bosonic degrees of freedom 
1731 [781, and the standard discrete-time quantum Monte Carlo 
action of Eq. ( [58] l is supplemented with the additional term 

GHiSOl: 


N Nr Nr 

/3^SB = -aEE E (61) 

2—1 T —1 sin 1 |t ) 

where a > 0 is the system-bath strength. This term intro¬ 
duces (non-local) ferromagnetic couplings between all spins 
along the Trotter direction of the dual classical spin system. If 
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FIG. 12. Simulation results for the final-time (t = f/) ratio Pi/Pc using SQA with explicit baths of varying system-bath interaction strengths 
a, as defined in Eq. For all sweep values, the ratio grows as a is increased, eventually crossing or saturating at Pi/Pc = 1, denoted by 
the dashed line. SQA simulation parameters: /3 = 10, Nr = 64, a linear annealing schedule (A(s) = 1 — s, P(s) = s), and a total of (a) 100, 
(b) 200, and (c) 500 sweeps. 10® SQA runs were performed for each a value. Error bars were generated by performing 100 bootstraps on the 
10 ® runs and taking twice the standard deviation of the Pi/Pc values. 


this coupling overwhelms the transverse field term in the ac¬ 
tion, which is the only other source of couplings in the Trotter 
direction for the transverse Ising Hamiltonian, then it is con¬ 
ceivable that all the spins in the imaginary-time direction be¬ 
have as one large spin flipping together. In this sense, updates 
in the dual classical spin system should behave as SA updates 
on the Ising part of the Hamiltonian. 

In order to test this intuition, we use the 8-qubit “quan¬ 
tum signature” Hamiltonian proposed in ref. ||24l, depicted 
in Fig. [TT] This Ising Hamiltonian has the feature of having 
a 17-fold degenerate ground state, 16 of which are connected 
via single spin flips, so they are called the “cluster” ground 
state. The remaining ground state is at least 4 spin flips away, 
so is referred to as the “isolated” ground state. We denote 
the average population in the 16 cluster ground states by Pc 
and the population in the isolated ground state by Pj, both 
at the final time t = tj. In the thermal equilibrium state 
all ground states are equally probable, i.e., Pi/Pc = 1. In 
Ref EM it was shown that SA and quantum annealing can be 
differentiated by the value of the ratio Pi/Pc'- SA will always 
have Pi/Pc > 1, i.e., SA preferentially populates the iso¬ 
lated ground state relative to any given cluster ground state, 
while quantum annealing will typically have Pi/Pc < 1, i e., 
quantum annealing preferentially populates the cluster ground 
states. These differences can be understood by studying the 
corresponding classical and quantum spectra. The classical 
spectrum is such that any random state can reach the ground 
states without encountering any local minima, and there are 
more paths to reach the isolated ground state than any given 
cluster ground state. This explains why SA favors the isolated 
ground state. Using first order perturbation theory, the degen¬ 
eracy of the ground states is broken by the introduction of the 
transverse field, and the isolated state has no overlap with the 
perturbed ground state. Hence, quantum annealing will genet¬ 
ically not populate the isolated ground state. More details can 
be found in Refs. Il24ll28l . 

If SQA with an explicit bath can interpolate between QA 
and SA by tuning the system-bath strength, it should be able 
to interpolate between the two regimes of Pi/Pc < 1 and 
Pi/Pc > 1 by increasing the system-bath coupling strength. 


Our simulations validate this picture, as shown in Fig. [T^ In¬ 
creasing the system-bath strength raises the Pi/Pc ratio from 
close to 0 to above 1, thus causing the final distribution of 
ground states to increasingly resemble the classical limit. The 
intermediate coupling strength regime is the one that cannot 
be captured by the WCL. 

We further test the dependence on the total evolution time, 
measured in terms of the the number of sweeps, where a 
sweep is one complete Monte Carlo update of all spins. In the 
case of a small number of sweeps [see Fig. 12(a)| , we observe 
the ratio first peaking and then decreasing towards 1. This 
shows, as expected, that for short evolution times the system 
requires a strong system-bath coupling in order to reach its 
equilibrium state (with Pi/Pc = 1), i-e., that the system equi¬ 
librates faster as the system-bath strength is increased. When 
the number of sweeps is high [see Figs. 12(b)| and |12(c)| , 
the non-monotonic behavior is replaced with a monotonic ap¬ 
proach towards the thermal equilibrium value of Pi/Pc = 1- 
In all cases shown here, the ground state probability is close 
to one. 

We can understand the WCL as being able to capture the 
Pi/Pc ^ 1 regime. Values of Pi/Pc away from this, as seen 
in Fig. [T^ represent regimes that cannot be captured by the 
WCL, and are thus of particular interest for SQA methods that 
include an explicit bath dependence, as done here. 


IX. CONCLUSIONS 

In this work we revisited the problem of decoherence in adi¬ 
abatic quantum computation (AQC), using a master equation 
approach and a quantum Monte Carlo approach. We argued 
that the common perception that decoherence is always detri¬ 
mental for quantum computation should be qualified in the 
adiabatic case, since the extent of decoherence-induced dam¬ 
age depends on whether the adiabatic quantum computer op¬ 
erates in the weak or singular coupling limit (WCL or SCL), 
or perhaps in an intermediate regime. A well-engineered de¬ 
vice should operate in the WCL, whence dephasing occurs in 
the system’s instantaneous energy eigenbasis. This is a form 
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of decoherence that preserves the coherence of the instanta¬ 
neous quantum ground state (and all other instantaneous en¬ 
ergy eigenstates). Therefore this form of decoherence, even if 
extremely rapid, does not necessarily negatively impact AQC. 
This is in stark contrast to dephasing in the circuit-paradigm 
of quantum computing, where dephasing in the computational 
basis spoils the efficiency of the computation and mandates er¬ 
ror correction. Of course, there is no free lunch in the AQC 
either, as the WCL does allow decoherence in the form of ther¬ 
mal excitations, which can depopulate the ground state and 
thus negatively impact AQC. However, the final state reached 
in the long-evolution time limit is the thermal Gibbs state, and 
as long as the ground state population does not decrease with 
increasing system size, AQC can succeed. Whether a quan¬ 
tum speedup is possible is a separate question that depends on 
how rapidly the Gibbs state is reached, and which we have not 
addressed in this work. 

A particularly noteworthy result we have demonstrated here 
(see Fig. is that there is an optimal (problem-dependent) 
evolution time t f that is much shorter than the adiabatic time- 
scale, where the ground state population can be significantly 
higher than that in the thermal Gibbs state. We have also pre¬ 
sented evidence that the optimal tf decreases with the strength 
of the system-bath coupling (see Fig.|^. This shows that it can 
be advantageous to run AQC with a much shorter duration 
total evolution than that suggested by the standard (heuris¬ 
tic) inverse gap criterion. Finding this optimal evolution time 
might be done, e.g., by picking some initial tf, finding the 
ground state probability, then repeating the experiment using 
both Q/2 and 2tf, etc., thus performing a grid search that will 
converge rapidly on the optimal tf. Another interesting pos¬ 
sibility is to apply the theory of optimal stopping times for 
continuous processes EH- 

It is also interesting to note that the standard strategy of sup¬ 
pressing detrimental thermal excitations by increasing the sys¬ 
tem’s energy gap is only guaranteed to work in the very large 
gap limit, since as we have demonstrated using an Ohmic bath 
model, the excitation rate first rises and only then decreases 
as a function of the energy gap. This is true also in the multi¬ 
qubit setting [Eq. ([56])], since the suppression of thermal ex¬ 
citations via the gap-dependent Boltzmann factor is counter¬ 
acted by the gap-dependent excitation rate 7. These conclu¬ 
sions apply under the assumption that the Markovian adiabatic 
master equation we have used here holds; the situation in the 
non-Markovian setting can be different, and even an infinite 


energy gap may not suffice ll82]| . 

In stark contrast to the WCL results, in the case of the 
SCL decoherence occurs in the computational basis, resulting 
in the loss of instantaneous quantum ground state coherence. 
AQC then becomes impossible since the resulting final state is 
essentially fully mixed, i.e., the ground state population drops 
exponentially with system size. 

The WCL and the SCL describe two dynamical limits of 
the system-bath coupling, and we have described an approach 
(SQA-EB) that applies in the intermediate coupling regime. 
This is accomplished by integrating out the bath degrees of 
freedom in the path integral formalism, resulting in an effec¬ 
tive action for the spin system, that can be used with quantum 
Monte Carlo methods to sample from the instantaneous Gibbs 
distribution. A method such as SQA-EB for simulating open 
system quantum annealing beyond the WCL and the SCL is 
particularly useful in studying the effect of decoherence, as it 
provides us with a mechanism to controllably decohere qubits 
into classical bits. Although SQA-EB does not capture unitary 
dynamics during the quantum evolution, it has the advantage 
of enabling the study of large system sizes. This provides 
us with a valuable tool for modeling large open quantum an¬ 
nealing systems. It is particularly important as a means to 
go beyond the WCL, whose validity depends on the system’s 
energy gap remaining large compared to the system-bath cou¬ 
pling strength [Eq. (|^], an assumption that will be violated 
for sufficiently large systems encoding computationally inter¬ 
esting ground states, whose gaps are known to shrink rapidly 
as a function of system size. Methods such as SQA-EB will 
play an important role in deciding the relative advantage of 
quantum annealing over classical annealing 1831 . An out¬ 
standing open problem is to develop a method that generalizes 
the WCL and SCL adiabatic master equations and can capture 
adiabatic quantum dynamics in the full intermediate coupling 
regime. 
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Appendix A: Off-diagonal components and decoherence in the 

WCL 


second term in Eq. ( |A3c| i. 

^^7q/3(0) Lpfi^acLl^ O cbPcc 
Oi(3 c^a,b 


(A3a) 


We would like to isolate the contribution of the diagonal el¬ 
ements of p to the off-diagonal term (a| £wcl[/3] \b) given in 
Eq. ( |53| l. The reason is that we would like to demonstrate that 
under the right assumptions (non-degeneracy, as explained be¬ 
low), the diagonal elements of p do not appear in this off- 
diagonal term, which means that the off-diagonal elements of 
p evolve independently from its diagonal elements. This will 
allow us to extract the T 2 time. 


In addition to the KMS condition 7q,/3 (— w) = e ^“7/3a (w) 
we repeatedly use the following easily verihed identities; 


^^7a/3(0) I o,ab 


1 Paa 

aP \ 


1 

(A3b) 

^^7q/ 3(0) I ^/3,0,abP^ o.bb ~ 

0 ^ ^ ^Q(,0,ac'^/3,0,cfe , 

I Pbb 

aP \ 

/ 

(A3c) 


We repeat the same process for w 7 ^ 0 and use the KMS con¬ 
dition, but now all the diagonal terms (with equal Roman sub¬ 
scripts on the L operators) vanish due to Eq. ( |Alb[ ). With this 
in mind Eq. (|A4[l can be read off directly from Eq. (|A3|l. 



(Ala) 

'a,CL»7^0,aa — 0 

(Alb) 

'a,CL),ab — u) ba — ^uj,eb — Sa^ct^ab 

(Ale) 


(Aid) 


Eq. ( |Ala| i follows from Eq. ([^ provided the operators 
are Hermitian, which can be assumed without loss of gen¬ 
erality. Eqs. ( |Alb[ ) and ( |Alc[ l follow directly by taking ma¬ 
trix elements of La^uj- Eq. ( |Ald| ) follow since the term 
ac is non-zero only if the energies associated with the 
states {|£a)) kc)} satisfy oj = £c — £a', similarly, the term 
^ is non-zero only if w = £c — £&. Therefore the term 

Lp^u),acL^a lAi cb i^ non-zero only if £„ = Sb, i.e., the states are 
degenerate. 


Eor convenience we reproduce Eq. (|5T| here in three parts; 




Cc'>0 Q;/3 


c^a,b 


cj>0 aj3 


cj>0 aj3 


c^a,b 


c^a,b 


T t 

,oj,ac-^ct,u},cb 



(A4a) 

+e 



,cb) 

Pcc 

(4... 

ac^P,u>,cb 


(A4b) 

+e~ 

(SCJ T 

rt 

uj,ac-^^a,iu 


Paa 

(4,., 

ac^P,d),cb 


(A4c) 

+e~ 

(3cj T 

rt 

uj,ac-^ct,io 

,cb j 

Pbb 


We now make the simplifying assumption that no two states 
are degenerate, i.e., £a 7 ^ £b Va 7 ^ h. It then follows from 
Eq. ( |Ald| i that all terms in Eqs. ( |A3| l and ( |A4| i vanish since 
they are composed entirely of products of operators of the type 
appearing in Eq. (|Ald|l. 


(a £wcl[p] \b) = Pi + P2 + P3 

(A2a) 

-Pi = E E E lapi^)Lp,ui,acPcdLa,di.db 

(A2b) 

UJ Q;,/3 c,d 

P2 — „ E! E! E! lo‘Pi'^)Pac^''a.ui.cd^P,^,db 

(A2c) 

OJ a,l3 c,d 

P3 — 2 ^ E E ]^0!p{^)Lti.ui.ac^P,ai,cdPdb 

(A2d) 


OJ q;,/ 5 c,d 


Considering Eq. ( |A2| l, there are three opportunities for the 
diagonal elements to appear, in pd with c = d (Pi), in 
with c = a (P 2 ), and in p^b with d = b (P 3 ). Thus we 
consider these separately. Setting c = d and w = 0 in Pi 

yields Eap Ta/ 3 ( 0 ) Ec ;^M,“'=-^«.o,cbPcc- We split this into 
c 7 ^ a,b, which gives ( |A3a| l; c = a, which gives the hrst 
term in Eq. ( |A3b[ ), and c = b, which gives the first term in 
Eq. ( |A3c| ). Setting c = a and w = 0 in P 2 gives the second 
term in Eq. (|A3b|). Setting d = b and a; = 0 in P 3 gives the 


Therefore, we have shown that in the absence of degen¬ 
eracies the equations for the dissipative dynamics of the off- 
diagonal density matrix elements do not involve diagonal ele¬ 
ments of the density matrix. We thus consider the contribution 
of the off-diagonal elements next. To this end we set c 7 ^ c? in 
Pi, c 7 ^ a in P 2 , and c? 7 ^ 5 in P 3 . This yields, respectively; 

'7q/3 ,u;,ac PcdL^a,ui,db (A5a) 

<jj cx.,0 c^d 

^^ = -2EEEE7«.( ^)Pac^(y^Qj^c(l^0,uj,db (A5b) 

uj a,/3 c^a d 

^3 = -2EEEE7 ap{^)pDi,u],ac^l3,‘^,cdPdb (A5c) 

uj a,P c d^b 

Considering P{, using Eq. ( |Alc| ) we have Lp,u,,acL^a,u,,db = 
5uj,ea-eJu,ei-eiAp^acAa,db- Likewise, considering P^ 

we have d^a,uj,cd^P,^,db — —— 

5u,et-eJet,eAa,cdAp^db, and Considering we have 

Py.,u),ac^P,^,cd — d^,ea—E:cbuj,€d — £cAa,acAp^cd — 
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6u,,eo.-eJs^,ea^a,acAp^cd- In the absence of degenera¬ 
cies, we can further simplify by using that = 6b,c and 

6^ = Sa d- This then gives simplified expressions for P[, 

p^Pk- 


Eq. ( |A 6 a| i is proportional to 7 ^^ (0). After extracting the terms 
proportional to 7 q/ 3 ( 0 ) from P 2 and Pg this then reduces the 
dissipative contribution to the equation for the off-diagonal 
density matrix elements to; 


Pi ~ ^ 'Japi^d ^b)Sec — ea,ed — £b^P,ac^a 

i,dbPcd 

a,p c^d 

(A 6 a) 

P2 — 2 ^ ^ ^ ^ 'ya.pi^^b ^d^-^a.,bd-^p,dbPab 

(A 6 b) 

a,/3 d 


P3 ~ 2 ^ ^ ^ ^ 'yapi^a ^c)-^a.,ac^p,caPab 

(A 6 c) 

a,p c 



(aI Ewcl[p] I^) = 


Q:,/3 


7o/3(0) ( A /3,aa-^a,bb 2 


2 I 2 ^ ^ 'ya^i^b ^a')^<^,ba' -^^,a'b 

^ a'^b 


cy ^ ^ 'TcK/3(‘^a ^b'^^a,ab' 


j3b' a 


b'^a 


Pab 


(A 8 a) 


We next consider the contribution of the unitary dynamics, 
—i (a| [Hs{t) + HLs{t): p{t)] \b) [including the Lamb shift, 
Eq. ([Uj]. By a similar argument to the one above, it is not 
hard to see that this too only involves off-diagonal density ma¬ 
trix elements. Then, after combining with (a| £wcl[p] \b) = 
P[ -f P 2 + Pg, we finally obtain for the off-diagonal elements: 


= -H^a -eb)Pab 
at 

^ ^ ^ ^ ^ (^g ^b')-^a,ab'-^p^b'aPab 

a,/3 _ b' 


+E 


^ ^ Spipif^b ^a'^-^ct^ba'-^p,a'bPab 
a' 

^ ^ 'yapi^^b' —eai^h' '^OL^ybPa'b' 


a'^b' 

~2 E/ 7a/3(£b 

a' 

-- lapi^a 
b' 


^a' ^A(^,ba' A^ d'bPab 
A^^b'aPab 


(A7) 


The important point is that only off-diagonal elements of the 
density matrix appear in this expression. By complete posi¬ 
tivity of the Lindblad form of the adiabatic master equation 
na, the instantaneous eigenvalues of this system of linear 
equations must have magnitude less than or equal to zero. It 
follows that the off-diagonal elements of the density matrix 
all decay to zero, and hence that coherence between instanta¬ 
neous energy eigenbasis states decays away. (The case of zero 
eigenvalues correspond to stationary states, which include the 
instantaneous Gibbs state and any decoherence-free subspace. 
Eor the case of independent decoherence, where each qubit in¬ 
teracts with a separate bath, the latter is not possible.) 

The system equations ( |A7| i is coupled due to the appear¬ 
ance of pa'b' in the fourth line, and this prevents us from eas- 

(e) 

ily extracting the T 2 time. This difficulty arises from the P[ 
term in Eq. ( |A 6 a| i. We can make analytic progress in find¬ 
ing the time if further assume that not only the energies 
are non-degenerate, but also the energy gaps. I.e., we assume 
that 5E^-ea,ed-eh ^ 6a,c6b,d- Under this assumption P{ in 


(e) 

The T 2 time can now be simply read off to give the expres¬ 
sion in Eq. 


Appendix B: Ground state population loss in the WCL 


Starting again from Eq. ( p2) , we can write; 
(a I £wcl[p] I a) = 


uj a,p 


'y ] ( Pp,‘^,acPcdPa,u],da r)Pa,u:,acPP,‘^,cdPda 

c,d ^ 


2PcLcLa,uj,cdPp>‘^,da ) . 


(Bl) 


Eollowing a procedure similar to that used in AppendixfA] we 
note this expression can be significantly simplified in the ab¬ 
sence of energy and gap degeneracies. Under this assumption, 
we have, combining Eqs. ( |Alc| l and ( |Ald[ ): 

Pp,<^,acPa,dj,da ~ Aa,caAjd,ac6u,£c — Sa^dc : (B2) 

P^a,uj,acPP,cd AEi,acA^,c,a6uj,ea—^c^o,d i (B3) 

Pa u cdPp^^^da — Aq; ,ad-^P,da^uj,6a—€d^ca ■ (B4) 

Using these identities, we can simplify Eq. •ET) to; 


I ^wcl[p] I a) = EE 


'yapi^c ,acAoi,caPc 


^ap(,^a ^c')Aa,(icAp,caPaa\ 

^ ^ ^ [7a/3(^c ^a)Ai^,(icAa,caPc 


c^a a,P 


^ap{_^a ^c)AQi^,acAjd,caPaa\ j 


(B5a) 


(B5b) 


where in the second equality we used the fact that the c = a 
terms cancel. The dissipative dynamics of the population of 
state a is now in the form of a rate equation, with the positive 
terms representing repopulation of the state a, while the nega¬ 
tive terms represent depopulation of the state a. Therefore we 
can identify the rate of loss of population from the ath energy 

















21 


eigenstate due to the dissipative dynamics as 

OL^{,^a ^c)'^a,ac-^/3 ,ca • (B6) 

c^a 

We obtain Eq. (|5^ upon restricting to the ground state (a = 0) 
and assuming independent baths [ 7 ct/ 3 (w) = (5a/37(w)]. 


